Accepted for publication in the Astrophysical Journal 
Preprint typeset using LAT^^ style cmulateapj v. 11/12/01 



STELLAR KINEMATICS AND METALLICITIES IN THE LEO I DWARF SPHEROIDAL GALAXY 

- WIDE FIELD IMPLICATIONS FOR GALACTIC EVOLUTION 

Andreas Koch^'\ Mark I. Wilkinson^, Jan T. Kleyna^, Gerard F. Gilmore^, 
Eva K. Grebel\ A. Dougal Mackey^, N. Wyn Evans^, and R. F. G. Wyse'' 

akoch@astro.ucla.edu 
Accepted for publication in the Astrophysical Journal 

ABSTRACT 

We present low-resolution spectroscopy of 120 red giants in the Galactic satellite dwarf spheroidal 
(dSph) Leo I, obtained with the GeminiN-GMOS and Keck-DEIMOS spectrographs. We find stars with 
velocities consistent with membership of Led out to 1.3 King tidal radii. By measuring accurate radial 
velocities with a median measurement error of 4.6 kms^^ we find a mean systemic velocity of 284.2 kms~^ 
with a global velocity dispersion of 9.9 km s^^. The dispersion profile is consistent with being flat out to 
the last data point. We show that a marginally-significant rise in the radial dispersion profile at a radius 
of 3' is not associated with any real localized kinematical substructure. Given its large distance from 
the Galaxy, tides are not likely to have affected the velocity dispersion, a statement we support from a 
quantitative kinematical analysis, as we observationally reject the occurrence of a significant apparent 
rotational signal or an asymmetric velocity distribution. Mass determinations adopting both isotropic 
stellar velocity dispersions and more general models yield a M/L ratio of 24, which is consistent with the 
presence of a significant dark halo with a mass of about 3 x 10^ Mq, in which the luminous component 
is embedded. This suggests that Leo I exhibits dark matter properties similar to those of other dSphs 
in the Local Group. Our data allowed us also to determine metallicities for 58 of the targets. We find 
a mildly metal poor mean of — l.Sldex and a full spread covering Idex. In contrast to the majority of 
dSphs, Leo I appears to show no radial gradient in its metallicities, which points to a negligible role of 
external influences in this galaxy's evolution. 

Subject headings: galaxies: dwarf — galaxies: individual (Leo I) — galaxies: kinematics and dynamics 

— galaxies: evolution — galaxies: abundances 



1. INTRODUCTION 

Leo I is the most remote dwarf spheroidal (dSph) galaxy 
generally believed to be associated with the Milky Way 
(MW) subsystem in the Local Group (see Table 1 in 
Grebel, Gallagher, & Harbeck 2003). Due to its remark- 
ably high radial velocity of about 287 km s~^ (Zaritsky 
et al. 1989; Mateo et al. 1998) and its current large dis- 
tance of ~ 254 kpc (Bellazzini et al. 2004), Leo I is a crucial 
tracer for determining the total mass of the Galaxy (e.g., 
Kochanek 1996). The question of whether Leo I is actually 
bound to the Galaxy remains open (e.g., Byrd et al. 1994). 
A quantitative assessment of the possible role of Galactic 
tides in the internal (dynamical, star-formation history, 
etc.) evolution of the Galactic satellite dSphs is of consid- 
erable importance, for its implications for the evolution of 
small galaxies, and for the distribution of dark matter on 
small scales. In particular, Leo I is of great interest here, 
given its very large Galactocentric distance, as it is the 
least-likely dSph to be affected by tides (but see also the 
recent work of Sohn et al. 2006, hereinafter S06). Simi- 
larities in its dynamics with those of more nearby dwarfs 
would strengthen our confidence that Galactic tides are 
not a major effect in determinations of the dark matter 
distribution both of Leo I itself and the other dSphs 

The multiplexing capabilities of present-day multi- 



object spectrographs on 4-10 m class telescopes have 
opened a window to study dSph galaxies in unprecedented 
detail by enabling one to gather data sets for a large num- 
ber of individual stars. As a result, the systemic velocities 
of the majority of the Local Group (LG) dSphs are known 
(tabulated in the reviews by Grebel 1997; Mateo 1998; 
van den Bergh 1999; and more recent measurements as 
quoted in, e.g., Evans et al. 2000, 2003). Moreover, high- 
precision radial velocity data for many individual stars in 
all the known dSphs are being obtained. These serve as 
the kinematical input for detailed dynamical studies (Ma- 
teo et al. 1998; Cote et al. 1999; Kleyna et al. 2001, 2003, 
2004; Lokas et al. 2002; Wilkinson et al. 2004; Tolstoy et 
al. 2004; Chapman et al. 2005; Muiioz et al. 2005, 2006a; 
Walker et al. 2006a, 2006b; Wilkinson et al. 2006a,b; and 
references therein), which are determining the spatial dis- 
tribution of dark matter on the smallest available scales. 

DSphs are the least massive and least luminous galax- 
ies known to exist. They are characterized by absolute 
y-band luminosities My ^ —14 mag, low surface bright- 
nesses of >22magarcsec~^ and H I masses of less than 
10^ Mq (Grebel et al. 2003 and references therein). They 
have high central stellar radial velocity dispersions, typ- 
ically of order lOkms"^ (Aaronson 1983; Mateo 1998), 
which, together with their globular-cluster-like low lumi- 
nosities (Gallagher & Wyse 1994; Mateo 1998) and their 
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large core size - hundreds of parsccs rather than the few 
parsecs of star clusters - demonstrate that these galax- 
ies are dominated by dark matter at all radii (Mateo et 
al. 1997; Wilkinson et al. 2004, 2006a). This deduction 
presumes the dSph are in a state of dynamical equilib- 
rium, a presumption supported by direct evidence that 
the dSph are not highly elongated down our lines of sight 
(e.g., Odenkirchen et al. 2001; Mackey & Gilmore 2003; 
Klessen et al. 2003). This asumption merits careful test, 
however, especially in the outer-most parts of the dSph, 
where Galactic tides might become significant. 

The inferred high mass-to-light ratios (M/L), obtained 
under the assumption of dynamical equilibrium, are as 
large as 500 {M/L)q or more, making dSphs the small- 
est observable objec;ts available to study the properties of 
dark matter halos. Available measurements suggest that 
the dark matter in dSphs has some common properties, 
including a minimum halo mass scale of ^ 4 x 10^ Mq, a 
small range in central (dark) mass densities, and a mini- 
mum length (core radius) scale of order 100 pc (Wilkinson 
et al. 2006a). If established, these are the first character- 
istic properties of dark matter. 

Cold dark matter (CDM) simulations of galaxy forma- 
tion predict the existence of a large number of sub-haloes 
surrounding larger galaxies, such as the Milky Way or 
M31 (e.g., Moore et al. 1999). The correctness of this 
prediction remains a subject of debate: the number of ob- 
served dSphs is much less, though several have been dis- 
covered very recently, indicating that the known sample is 
incomplete (Belokurov 2006a, 2006b; Martin et al. 2006; 
Zucker et al. 2006a, 2006b, 2006c). The correspondence 
between prediction and observation more importantly re- 
quires masses, so that dynamical studies of the dSphs are 
critical. 

The spatial distribution of dark matter on small scales 

is an additional test of the nature of dark matter. Nu- 
merical simulations of dark matter halos in which it is as- 
sumed that there is no physical process except gravity act- 
ing remain limited in spatial resolution, but in general are 
characterized by a central density cusp p{r) oc r~° with 
a=l 1.5 (e.g., Navarro, Frcnk & White 1995, hereafter 
NFW; Ghighna et al. 2000). It is far from clear that these 
expectations are consistent with observations of galaxies 
significantly more massive than the dSphs (de Blok et al. 
2001; Spekkens et al. 2005), which are better accounted for 
by flat density cores. The available data on dSph kinemat- 
ics appear to be consistent with cored halos as well (Lokas 
2002; Read k Gilmore 2005; Wilkinson et al. 2006a). In 
the Ursa Minor dSph a central, constant-density core has 
been detected (Klcyna et al. 2003), and also in the Fornax 
dSph a core may be more likely than a cusped distribu- 
tion, although these inferences are still subject to large 
uncertainties (Strigari et al. 2006; Goerdt et al. 2006). It 
has been suggested that the cusp/core distinction could 
be the first indications of the properties of the physical 
particles which make up CDM (Read & Gilmore 2005), 
though it is essential that one allows for dynamical evo- 
lution associated with the astrophysical evolution of the 
dSph progenitors (Grcbel et al. 2003). 

According to a second scenario, the large observed ve- 
locity dispersions may be explained if dSphs are unbound 
tidal remnants that underwent tidal disruption while orbit- 



ing the Milky Way (Kuhn 1993; Kroupa 1997; Klessen & 
Kroupa 1998; Klessen & Zhao 2002; Fleck & Kuhn 2003; 
Majewski et al. 2005), and all happen to be aligned to 
mimic a bound system. Predictions following from this 
scenario include a significant depth extent along the line 
of sight (Klessen & Kroupa 1998; Klessen & Zhao 2002). 
However, in two of the closest Milky Way companions, the 
Draco dSph (Klessen et al. 2003) and the Fornax dwarf 
(Mackey & Gilmore 2003) direct evidence against such ef- 
fects is available. Other arguments against dSphs being 
unbound tidal remnants without dark matter include that 
they all experienced extended star formation histories with 
considerable chemical enrichment and that they follow a 
metallicity-luminosity relation (sec Grebel et al. 2003 for 
details). In the light of its large Galactocentric distance, 
Leo I is well suited as a test of the possible importance of 
tidal perturbation. 

The subject of tides in the context of dSph galaxies suf- 
fers from a very distorted literature. In large part this 
follows from an historical labelling issue. Specifically, the 
projected surface brightness distributions of dSph galax- 
ies are nowadays preferentialy fit with King models. A 
King model has two named linear scales: an inner scale, 
frequently called a 'core radius' and an outer truncation 
radius, called a 'tidal radius'. While these labels have 
some relevance in their historical application to models 
of globular clusters of stars truncated by Galactic tides, 
they have no general physical interpretation in a differ- 
ent context. In general, dSph stellar populations are not 
single-component, and there is no astrophysical reason to 
assume ab initio that the distribution function is simply 
isothermal. Some form of numerically tractable functional 
fit to a projected surface brightness profile is required in 
the dynamical analysis, but a more sophisticated analy- 
sis is required than fitting functions to surface brightness 
profiles before the possible importance of tides can be in- 
vestigated. It is sometimes suggested that high-frequency 
distorted photometric structure is indicative of tidal dam- 
age, but no available numerical simulations support this 
view: rather, if tides are important, they produce smooth 
systematic distortions in the outer parts of a system. One 
test, applicable in some circumstances, is detection of ap- 
parent rotation or spatially asymmetric velocity distribu- 
tions associated with such a distortion, as claimed to be 
present in Leo I by S06. We will consider that in our anal- 
ysis of Leo I. Some of these relevant issues are discussed 
further by Read & Gilmore (2005), and by Read et al. 
(2006a, 2006b). 

Based on its projected position on the sky, Leo I was 
suggested to be part of a Galactic Fornax-Leo-Sculptor 
Stream (Lynden-BeU 1982; Majewski 1994), i.e., that it 
lies along a Galactic polar great plane together with these 
other dSph galaxies. However, despite the lack of proper 
motion and orbital measurements for Leo I, the radial ve- 
locity and kinematical data to hand apparently rule out a 
physical association of Leo I with such a physical group 
(Lynden-BeU & Lynden-BeU 1995; Piatek et al. 2002, 
2006). The remote Leo I dSph holds a special place in 
the kinematic studies of MW satellites (Mateo et al. 1998, 
hereinafter M98). From detailed kinematical modeling of 
the velocity dispersion profile of 33 stars within the core 
region of Leo I, M98 concluded that this dSph contains a 
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significant amount of dark matter. 

Apart from its stellar kinematics, Leo I is interesting in 
terms of its star formation history (SFH; see also Bosler 
et al. 2006). The dSphs of the Local Group exhibit a 
variety of SFHs, where no two dwarfs are alike (Grebel 
1997; Mateo 1998). In spite of this diversity, all nearby 
galaxies studied in sufficient detail have been shown to 
share a common epoch of ancient star formation, although 
the fraction of old populations with ages > 10 Gyr varies 
from galaxy to galaxy (Grebel 2000; Grebel & Gallagher 
2004). Indeed Leo I has a prominent old stellar popula- 
tion (Held et al. 2000) , but more striking is the prevalence 
of intermediate-age populations - about 87% of its stars 
formed between 1 and 7 Gyr ago (Gallart et al. 1999a). 
The details of the SFH of Leo I were derived photometri- 
cally, but particularly in galaxies with mixed populations 
the age-metallicity degeneracy cannot be broken without 
spectroscopy. Our current inferences of the chemical evo- 
lutionary histories of nearby dwarf galaxies are still mainly 
based on photometry, while detailed and numerous spec- 
troscopy of as many stars as possible is invaluable for 
helping to securely constrain the SFHs or to quantify age- 
metallicity relations. It is also essential to study the chem- 
ical evolution of dSphs in order to understand their com- 
plex SFHs in terms of regulating processes such as the 
inflow or accretion as well as the outflow of gas (Carigi et 
al. 2002; Dong et al. 2003; Lanfranchi & Matteucci 2004; 
Hensler, Theis, & Gallagher 2004; Robertson et al. 2005; 
Font et al. 2006). We are currently conducting a program 
to measure the kinematics and chemical abundances in 
several Galactic dSphs and to correlate these with their 
evolutionary histories (e.g., Wilkinson et al. 2006b, Koch 
et al. 2006a). 

Here we present a kinematic and abundance study of 
the distant Galactic dSph Leo I. In addition to a detailed 
kinematic analysis, appropriate to the size of the available 
dataset, we have measured the average metallicity and the 
spread and shape of Leo I's metallicity distribution func- 
tion (MDF), which provide valuable insights on any possi- 
ble environmental dependence of the processes governing 
the evolution of dSphs. This paper is organized as follows: 
§2 describes our data acquisition and reduction and the 
determination of stellar velocities. In §3 we derive Leo I's 
global velocity distribution from which we investigate po- 
tential velocity gradients in §4. Section 5 is then dedicated 
to the analysis of Leo I's radial velocity dispersion profile 
and the influence of binaries on this profile. This radial 
profile allowed us to obtain mass and density estimates of 
this galaxy in §6. The role of an anisotropic velocity dis- 
tribution on dispersion and mass estimates is addressed in 
§7. In §8 we derive metallicities for our targets and dis- 
cuss the implications for Leo I's chemical evolution, and §9 
finally summarizes our results. 

2. DATA 

2.1. Target selection 

Targets were selected from photometry obtained within 
the framework of the Cambridge Astronomical Survey 
Unit (CASU) at the 2.5 m Isaac Newton Telescope (INT) 
on La Palma, Spaii|3. From this, we selected red giants so 
as to cover magnitudes from the tip of the red giant branch 

^ see |http : //www ■ ast ■ cam. ac .uk7$\sim$wf csur/^ 



(RGB) at / 18 mag (Bellazzini et al. 2004) down to 
~1.6mag below the RGB tip, going as faint as / < 19.6. 
At our faintest observed magnitudes, our spectra reach 
signal-to-noise (S/N) ratios of ~5, which is sufficient to 
enable accurate radial velocity measurements at our low 
resolution. 

2.2. Observations 

One part of our observations was carried out in queue 
mode with the Gemini Multiobject Spectrograph GMOS 
at the 8-meter Gemini North telescope, Hawaii, over 
four nights between December 2004 and January 2005. 
These data were obtained under photometric conditions 
and a clear sky. The camera consists of three adjacent 
2048x4068 CCDs, which are separated by gaps with a 
width of ~37 pixel each. We used 0'.'5 wide multiob- 
ject slits and the R400-I-G5305 grating set in two modes, 
with central wavelengths of 8550 A and 8600 A, provid- 
ing a spectral coverage of 1800 A. This dithering enabled 
us to efficiently interpolate the spectra across the CCDs' 
gaps during the later coaddition process. The CCDs were 
binned by 2x4 during readout and the resulting spectra 
have a nominal resolution of 2200. In total, we targeted 
68 red giants spread over three different fields in the east- 
ern half of the galaxy (see Fig. 1, Table 1). Each of the 
three fields was exposed for three hours in total, where we 
split up the exposures into 3 x 1800s per setup (i.e., each 
of the two wavelength modes) to facilitate sky subtraction 
and cosmic ray removal. Due to bad weather, the central 
pointing (Field 3, see Table 1) was exposed for only 3600 s 
in total. 

The other part of our data was obtained with the 
DEep Imaging Multi-Object Spectrograph (DEIMOS) at 
the 10-meter Keck II telescope. Leo I was observed on 
the DEIMOS multislit spectrograph on the Keck II tele- 
scope on February 6 to February 8, 2005. Conditions were 
mediocre, with high humidity and poor seeing, but ap- 
proximately six hours of Leo I observations were obtained, 
consisting of five slit masks concentrated on the dSph's 
periphery. We used the 1200 line/mm grating. The slits 
were 0.7" in size, and our resolution was about 0.3 A per 
pixel. 

2.3. Data Reduction 

2.3.1. GMOS data 

The GMOS data were reduced using the standard 
GMOS reduction package in IRAF, starting with the bias 
frames taken as part of the day-time calibration process. 
Spectroscopic fiat-field exposures were obtained adjacent 
to each set of science exposures, from which we could 
model and remove the spectral signature and spatial pro- 
file of the illumination pattern. Through dividing by the 
remaining normalized flat-field frame, we ensured a com- 
plete removal of any residual sensitivity variation. Since 
the CuAr wavelength calibration frames were taken up to 
four days before and/or after the science spectra, during 
which time the telescope flexure may significantly change, 
we had to rely on the skylines of each individual science ex- 
posure to obtain an accurate wavelength solution. Hence, 
we chose a number of strong and isolated, unblended night- 
sky OH-emission lines around the CaT region from the 
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atlas of Osterbrock et al. (1996, 1997). Fitting the re- 
spective lines with a low- order polynomial yielded a RMS 
uncertainty in the wavelength calibration of the order of 
0.05 A, corresponding to ~ 2kms~^ at our region of in- 
teresl[3. A crucial step in the analysis of CaT absorption 
features consists of the accurate subtraction of the promi- 
nent skylines in the near-infrared spectral region. While 
the centers of the first two of the Ca triplet lines are un- 
affected by adjacent skylines at Leol's systemic velocity, 
the third Ca absorption line coincides with the dominant 
OH-band at ~8670A. These undesired contaminants are 
efficiently removed by the software by fitting a low order 
polynomial to the background and subtracting this fit col- 
umn by column from the science spectra. Moreover, the 
accuracy of our sky subtracted spectra (after coaddition) 
was improved by the fact that we obtained spectra with 
the slit centred both at the centre and offset towards the 
right and left hand side of the stellar seeing disk, yield- 
ing a better coverage of the average sky background. Fi- 
nally, the extracted spectra were shifted towards the local 
standard of rest barycenter and median-combined using 
a standard cr-clipping algorithm. The average S/N ratio 
achieved throughout our reductions is ^28, with a mini- 
mum of 5.3 at I = 19.6. 

2.3.2. DEIMOS data 

The DEIMOS data were reduced using the DEEFQ 
pipeline, developed by the DEEP2 project at the Univer- 
sity of California-Berkeley and based on the IDL spec- 
tral reduction pipeline of the Sloan Digital Sky Survey 
(Abazajian et al. 2004). The DEEP pipeline first identifies 
spectra using tungsten lamp fiatfield exposures, and then 
finds a two dimensional (pixel-by-pixel) wavelength solu- 
tion for each spectrum, assisted by a DEIMOS-specific op- 
tical model. Typically, the solution is accurate to 0.005A. 
All calibration images were taken at the start of each night; 
this is possible because the spectrograph optical system 
is stabilized using an active flexure compensation system. 
Sky subtraction is performed using a 6-spline sky model. 
In the last stage, the sky-subtracted two-dimensional spec- 
tra are extracted into one-dimensional flux and variance 
spectra in the form of FITS tables. The variance spec- 
trum is flnally used to compute the spectral S/N ratio. 
The S/N achieved throughout our processing reaches ^-^45 
for the brighest targets and is as low as 5 for the faintest 
stars at I = 21.2. 

Generally, the DEEP pipeline worked as expected. For 
certain exposures, however, the slit identification and trac- 
ing in the fiats failed, and we found that the sensitivity 
parameter controlling slit detection had to be adjusted. 

2.4. Radial Velocities 

Radial velocities were derived from our final reduced set 
of spectra by cross-correlating the three strong Ca lines at 
AA8498, 8542, 8662 against a synthetic template spectrum 
of the CaT region using IRAF's cross-correlation package 
FXCOR. The template was synthesized using representa- 
tive equivalent widths of the CaT in red giants. The final 

^ Using the entire spectral region at hand, covering ca. 1800 A, yielded highly inaccurate results in the wavelength solution, as too large a range 
of distortions was being fitted with too few lines, resulting in off'sets between individual exposures of the order of ~ 50 — lOOkms"^. However, 
when concentrating on the narrow window (~500A) around the CaT, in which we are primarily interested, we circumvent any distortions 
towards the edges of the spectra and could thus obtain a highly improved accuracy. 

^ see http://astro.berkel6y.edu/~coop6r/d6ep/spec2d/ 



velocity difference between object spectrum and the rest- 
frame template was then determined from a parabolic fit 
to the strongest correlation peak. 

2.5. Velocity errors 

The formal velocity errors returned from the cross- 
correlation in IRAF have a median value of 1.7kms~^. 
However, we have found in past studies that the FXCOR es- 
timates are generally too optimistic (see e.g., M98; Kleyna 
et al. 2002). Since we have only single epoch data at hand, 
lacking the possibility to compare the velocity stability 
over a longer time period, we cannot estimate the true 
uncertainties through the discrepancies of measurements 
taken at different times (Vogt el al. 1995, M98). For the 
CMOS data, we rather follow the prescription of Kleyna 
et al. (2002) and divide the data of the brightest stars 
(with / < 19 mag) into two final sets of spectra, the first 
comprising those individual spectra with a central wave- 
length of 8550 A and the second being those centered at 
8600 A. Each of these sets was separately cross-correlated 
against the template. We use the discrepancy between the 
two resultant velocities as an estimate of the true measure- 
ment error. Hence, we rescaled the Tonry-Davis i?-based 
errors (Tonry & Davis 1979) returned by FXCOR to ob- 
tain the expected reduced discrepancy of unity. We 
found that we had to apply a scale factor of 2.7 to yield 
an accurate estimate of our velocity errors. In order to 
check whether our wavelength calibration, which is based 
on only a narrow spectral window around the CaT has 
introduced any bias, we also determined radial velocities 
by cross-correlating each of the three Ca lines separately 
against the same template. In this case, the mean devia- 
tion of the resultant velocities did not exceed 1.8kms~^. 
These uncertainties were added in quadrature to the for- 
mal measurement errors (including the above Tonry-Davis 
R scaling) to yield the final velocity error, which we show 
in Table 2. The final CMOS data set contained 68 stars 
with a median velocity error of 5.0kms~^. 

For the DEIMOS data, the nominal IRAF velocity er- 
rors were also rescaled by a common factor. The rescaling 
was performed by computing separate velocities for lines 
1+2 and line 3 of the CaT, and then requiring the differ- 
ence between the two velocities divided by the quadrature 
sum of the rescaled nominal errors to obey the expected 
statistics. The final DEIMOS member list includes 27 
member velocities with per pixel S/N > 10, and a median 
velocity error of 2.7 kms~^. 

3. RADIAL VELOCITY DISTRIBUTIONS 

To determine an overall mean velocity and a global dis- 
persion for Leo I, a prerequisite for an assessment of the 
targets' membership, we used the maximization method 
described in Kleyna et al. (2002) and Walker et al. 
(2006a). In this approach, the probability that an ob- 
served distribution of velocities Vi and corresponding mea- 
surement errors cr,- is drawn from a Gaussian distribution 
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with mean {v) and dispersion a is given by 



p{{vi,(Ti}\{{v),a)) cxn^i 



exp 



(1) 

It is mathematically more convenient to calculate the nat- 
ural logarithm of this expression, and maximizing Inp is 
equivalent to maximizing p itself. Rejecting stars deviat- 
ing more by than 3<t as likely non-members, we iterated 
until convergence (requiring that the solution from subse- 
quent iterations changed by no more than 0.005kms~^). 
The respective uncertainties in the mean and dispersion 
were then calculated from the covariance matrix (Walker 
et al. 2006a). This procedure yielded a mean systemic ve- 
locity of Leo I of (284. 2± 1.0) kms"-'^ and a dispersion of 
9.9±1.5kms-i from the combined CMOS plus DEIMOS 
sample, (284.6±1.5, 9.9±2.1) kms^^ for the CMOS data 
alone and (283.8±1.5, 9.9±2.1)kms~^ as derived from the 
DEIMOS targets only. The mean velocities and global dis- 
persions are tabulated in Table 3 separately for the differ- 
ent fields and Fig. 2 shows the resulting histograms. 

This compares to a value of 287.0 ± 1.9 km s""'^ as found 
by M98 and their respective dispersion of 8.8 ± 1.3kms~^ 
for stars in the innermost 3. '5. 

If we take the Sc-eut as a membership criterion, 5 of 
the 68 GMOS targets and 5 of the red giants from the 
DEIMOS sample have to be rejected as foreground con- 
tamination, where 4 out of these 5 stars in each sample 
have velocities below 200kms~^. In order to generate the 
final sample for determining the dispersion profile of Leo I 
we combined both of our GMOS and DEIMOS sets by 
shifting them to a common median velocity. This proce- 
dure is justified by the lack of any significant radial gradi- 
ent in the velocities from DEIMOS, which is found to be 
(—0.04 ± 0.18) kms~^ arcmin"^ so that such a shift does 
not introduce any artificial velocity gradient or a falsifica- 
tion of the dispersion profiles. At the end of this procedure, 
we found that our median velocity from the GMOS plus 
DEIMOS set differed from that of M98 by -2.8 kms"!, 
a difference significant only at the p = 0.31 level (note 
that our sample has no stars in common with the M98 
sample, which would allow for a direct comparison of the 
velocities). These values are finally listed in Table 2. 

4. KINEMATIC TESTS FOR TIDAL DAMAGE AND 
APPARENT ROTATION 

Those dSph galaxies studied so far show no significant 
contribution from angular momentum to the systems' dy- 
namical structure: dSphs obtain their size and shape from 
anisotropic random pressure. This is, incidentally, an ar- 
gument against an evolutionary transformation from dwarf 
irregular to dSph galaxies via ram pressure stripping (e.g., 
van den Bergh 1999; Grebel et al. 2003). More fundamen- 
tally, apparent rotation (generally visible in projection) is 
a characteristic signature of tidal disturbance. In this con- 
text, apparent rotation of the stellar component of a dSph 
is a test of possible tidal 'heating' of the stellar orbits in 
the course of the galaxy's orbit within the external tidal 
field of the Milky Way. N-body simulations, e.g., of Oh 
et al. (1995) have shown that the velocity dispersion of 
a tidally-limited system is sustained at the actual virial 
equillibrium value. On the other hand, these simulations 
suggest that any strong dynamical effects on a dSph from 



the tidal field of the Galaxy can invoke streaming motions 
^y. — ^^^)2i]« the outskirts of the dwarf, which will lead to a system- 
atic change in the mean velocity along its major axis, thus 



imicking rotation (Piatek & Pryor 1995; Oh et al. 1995; 
Johnston et al. 1995; M98). A more detailed analysis of 
tidal effects on dSph galaxies is provided by Read et al. 
(2006b). 

It has frequently been suggested in the literature that 
the detection of member stars well beyond the formal tidal- 
limit radius of dSphs, c;anonic;ally defined via a single- 
component King model representation of the observed sur- 
face brightness profiles, is evidence for physical tides. A 
large number of dSphs have been claimed to show hints 
of extratidal stars and thus tidal disruption (Irwin & 
Hatzidimitriou 1995; Martmez-Delgado et al. 2001; Palma 
et al. 2003; Muiioz et al. 2005, 2006a, 2006bj. Evi- 
dence for this has also been reported for Leo I (S06). What 
such studies show, in fact, is that the presumed paramet- 
ric fit to the measured surface brightness (and/or direct 
star counts) is inappropriate at large radii. Given that 
there is no astrophysical basis to the application of a King 
model to a galaxy, an astrophysical interpretation of the 
corresponding labels as physical processes is fraught with 
peril. Evidence that a specific dSph is being affected by 
tides suflaciently strongly that its kinematics are affected 
can be provided most efficiently from a kinematic analysis. 
A test for apparent outer rotation is such a test. 

With its present-day Galactocentric distance of 
(254±19)kpc and its high systemic velocity, Leo I is not 
expected to be affected by tides. M98 concluded from their 
kinematical data that the core of this galaxy has not been 
significantly heated by tides. S06 report on a large popula- 
tion of photometrical and radial velocity member red giant 
stars out to nearly two nominal tidal radii, which they in- 
terpret as evidence of tidal disruption. We note that also 
our sample comprises nine red giants outside rud, whose 
velocities suggest membership. We also emphasise that 
these statements tell us that the outer parts of Leo I are 
not well-fit by a single King model. Whether they tell us 
anything about the dynamics of Leo I remains question- 
able for the moment and has to await further dynamical 
tests. 

In order to assess whether there is any significant indi- 
cation of rotation in the Leo I velocities, which may help 
to elucidate the role of tides, we pursued two tests. First, 
we show in Fig. 3 the spatial distribution of our measured 
radial velocities (including M98's 33 data points) along 
their projected major axis radial distance. In addition, we 
calculated the mean velocity as a function of radius, which 
is shown in the lower panel of Fig. 3. It is worth notic- 
ing that this mean velocity remains fairly constant across 
the range covered by our data and the deviations from the 
global mean of our sample are well within the errorbars. 
Moreover, we do not detect any significant sign of a veloc- 
ity gradient in the inner regions as suggested by S06. If 
Leo I has been significantly affected by Galactic tides, one 
would expect an asymmetric excess of velocity outliers, 
with stars of higher or lower velocity typically lying on 
opposite sides of the galaxy due to the presence of possi- 
ble leading and trailing tidal tails. In this vein, S06 report 
on the presence of a number of stars with higher velocities 
than the average in the outer regions to the western half 
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of Leo I. The tidal model then predicts a comparable pres- 
ence of low velocity stars towards the eastern extension. 
Although our sample contains one red giant at a large 
(western) distance with a radial velocity that is higher by 
3.8 cT than the sample mean, there is no apparent tendency 
of a systematically lower mean velocity of the targets at 
large distances towards the east. 

While we find a mean radial velocity of 294.3±2.9kms~^ 
for the stars outside the formal tidal radius on the east- 
ern side, the value for the stars towards the western half 
including (rejecting) the high velocity outlier amounts to 
289.5±4.0 (285.7±5.5) kms'^ Prom this point of view 
our data do not show any hint of significant rotation, nor 
an asymmetry of the velocities on either side along the 
galaxy's major axis. 

Secondly, we calculated the difference in mean veloc- 
ity on either side of bisecting lines passing through each 
of the targets, assuming that each such line is a poten- 
tial rotation axis. The axis yielding the maximum mean 
velocity difference is then taken as the "rotation axis", 
with the associated velocity difference being an estimate 
of the amplitude of the rotation. Fig. 4 shows the result- 
ing velocity differences for our dataset. Due to the con- 
finement of our GMOS fields to the eastern half of Leo I, 
we cannot perform any test for rotation in this subsample 
alone. However, based on only the 56 candidate members 
from DEIMOS, we find a maximum velocity difference of 
1.2kms^^. Nevertheless, the spatial distribution of these 
stars does not allow one to unambiguously determine a po- 
tential rotation axis (see top panel of Pig. 4). In the case 
of the combined GMOS and DEIMOS sample, the ampli- 
tude of the potential rotation reaches a value of 5.2 km s"""^ 
(middle panel of Fig. 4). When also the 33 stars from 
M98 are added, the mean velocity difference changes to 
3.0kms~^ (bottom panel of Pig. 4). It is worth noticing 
that in both of the latter cases the axis of the maximum 
rotation signal occurs at a position angle (78°, 80°) which 
coincides with LeoFs position angle of (79 ± 3)° (Irwin & 
Hatzidimitriou 1995). 

In order to test the significance of these results, we con- 
structed 10'* random samples of radial velocities at the ob- 
served spatial positions. The velocities were taken from a 
normal distribution with the same mean and standard de- 
viation as in our observations, additionally allowing for a 
variation due to the measurement uncertainty. By means 
of this Monte Carlo simulation we find that the signifi- 
cance of the rotation in the DEIMOS sample is consistent 
with zero, whereas the amplitudes of the combined GMOS 
plus DEIMOS data sets and the full data including M98's 
stars amount to a significance of 86.7% (1.5 a) and 20% 
(0.3 (j), respectively. Thus the rotational signature we de- 
tected appears not to be statistically significant in either 
case. Consequently, we have no evidence of any kinemati- 
cal effect due to Galactic tides operating on Leo I. 

Another cause of apparent rotation is relative transverse 
motion. Unfortunately, there is no estimate of Leo I's 
proper motion available in the literature at present. Hence, 
we cannot rule out the possibility that the amount of ob- 
served rotation is entirely due to the relative motion be- 
tween the Sun and the dSph, Vrei{l, b), in the Galactocen- 
tric rest frame. A non-rotating object will appear to rotate 
as seen in the heliocentric rest frame, provided there is a 



significant gradient in this relative motion across the dSph 
(see Walker et al. 2006a). We can attempt the inverse cal- 
culation: we tested whether one can constrain Leo I's true 
proper motion by demanding that this galaxy does not 
rotate but the entire amplitude of the observed velocity 
gradient is caused by such projection effects. For this pur- 
pose, we calculated the Galactocentric rest frame velocity 
of each star through the formalism provided in Piatek et 
al. (2002) and Walker et al. (2006a). This was done for 
a wide grid of assumed proper motions, corresponding to 
transverse velocities from —500 to -1-500 km s~* in both ni 
and /Zfe. Regrettably, the full range of input {^i, m,) was 
able to reproduce the observed apparent rotation ampli- 
tude adequately, so no deductions are possible. We then 
re-ran the same formalism by adopting a toy distance of 
50kpc to Leo I and secondly assumed its line-of-sight ve- 
locity as zero. In both cases the parameter space of possi- 
ble proper motions increased even more. In consequence, 
we can neither unequivocally exclude the possibility that 
Leo I's observed, though insignificant, velocity gradient is 
purely caused by the aforementioned projection effects or 
due to its large distance or systemic velocity, nor can we 
conclusively restrict its proper motion to a reasonable es- 
timate. 

5. VELOCITY DISPERSION PROFILE 

From their sample of 33 stars, M98 reported cen- 
tral velocity dispersions ranging from (8.6±1.2) kms^* 
to (9.2±1.6)kms-i, which were obtained using several 
model-independent estimators. This is slightly lower than 
our global estimate from eq. 1, but consistent within the 
uncertainties. 

We determined the velocity dispersion profile using the 
maximum likelihood method outlined in Kleyna et al. 
(2004). The data were binned such as to maintain a con- 
stant number of stars per bin, where we chose the binsize 
such that no fewer than ten stars were included. A Gaus- 
sian velocity distribution around the single mean veloc- 
ity of the entire ensemble is then assumed for each radial 
bin, convolved with the observational errors. The mem- 
ber velocity distribution centred on the systemic velocity 
is then combined with an interloper distribution, Pini{v), 
contributing a fraction to each bin. This allows the 
probability distribution of the true dispersion a in the bin 
through an extension of eq. 1: 
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We then perform a maximum likelihood fit over a and 
fint^ marginalise over the interloper fraction and finally 
find the most likely dispersion from the marginalised dis- 
tribution. 

We take account of the effects of the Galactic foreground 

interlopers in our data set by assuming various associ- 
ated distributions Pint{v) over the whole range of interest 
around the systemic velocity of Leo I. Consequently, we 
tested power-law velocity distributions with varying in- 
dices against the uniform case and against each other and 
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found from a Kolmogorov-Smirnov (K-S) test that the re- 
sulting dispersion profiles were practically identical at a 
100% confidence level. The velocity distribution of the ten 
non-member stars in our sample itself is best represented 
by a power-law with an exponent of —0.7 (with a K-S prob- 
ability of 93%), whereas the K-S probability that these 
stars' velocities are uniformly distributed is 52%. We em- 
phasise, however, that the resulting dispersion profiles do 
not change with respect to the adopted parameterisation 
of the underlying interloper distribution. This is also to 
be expected, given the low total number of non-members 
sampled, a consequence of the overall high systemic veloc- 
ity of Leo I relative to line-of-sight galactic halo stars. We 
will, for convenience, use the uniform case in the following. 
Typically, the peak interloper fraction for the outermost 
bin does not exceed 0.1 and is compatible with zero in the 
inner bins. 

Finally, formal error bounds were determined by numer- 
ically integrating the total probability of the data set and 
finding the corresponding 68 per cent confidence intervals. 
Fig. 5 shows the resulting velocity dispersion profiles as 
a function of the projected spherical radius for different 
membership criteria, and separately for the subsamples. 

5.1. Radial variations 

Given the derived uncertainties in the profiles, the dis- 
persion profile of Leol can be considered to be approxi- 
mately flat out to the nominal tidal radius for the ma- 
jority of the data combinations shown. Exclusion of the 
3c7-outlier in the last bin leads to only a marginal change. 
Inclusion of the data from M98 with our own data into the 
maximum sample decreases the dispersion in the central 3' 
by ~2kms~^, while the associated rebinning of the data 
introduces a small drop in the dispersion beyond 10', quan- 
tifying the limited sensitivity of the binned dispersions to 
sample binning. 

5.2. Absence or presence of kinematical substructure 

Localized, kinematical structures in Local Group dSphs 

have been detected, specifically in the form of low veloc- 
ity dispersion (cold) spatially clumpy substructure in the 
UMi dSph (Kleyna et al. 2003; Wilkinson et al. 2004). In 
Sextans a low velocity dispersion (cold) central core has 
been found (Kleyna et al. 2004) and independently the 
presence of an off-centered kinematically distinct subpop- 
ulation has been reported (Walker et al. 2006b). Such 
distinct cold features can be interpreted as the remains of 
dispersing star clusters, with their central locations arising 
due to the clusters being dragged towards the inner regions 
by dynamical friction (Kleyna et al. 2003, 2004; Goerdt 
et al. 2006). An alternative proposal involves the after 
effects of a merger between a dSph and an even smaller 
system (Coleman et al. 2004). 

Despite the overall flatness of Leol's velocity disper- 
sion towards larger radii, the profile does show some radial 
structure. A noteworthy radial feature is an apparent rise 
in the dispersion profile at ~ 3' (r/rud ~ 0.25). This trend 
is persistent, regardless of radial binning and the chosen 
sample, and it is also inherent in the profile observed by 
Mateo (2005). The major difference between our profile 
and the one reported in Mateo (2005) is that his disper- 
sion rises continuously out to 10.'5 after having reached its 



minimum after the bump. Hence, it cannot be excluded a 
priori that the local dispersion max;imum may in fact be 
localized and physically real. 

In order to test whether our data set supports the pres- 
ence of any underlying localized kinematical substructure, 
we estimate the position-dependent dispersion using a non- 
parametric approach. We define the locally weighted, av- 
erage dispersion as 

Z^i=i K{x,Xi,y,yi,h) 
(Walker et al. 2006b, after Nadaraya 1964; Watson 1964). 
Here, Vi and (Tj are as usual the observed radial velocities 
and associated uncertainties of the N targets at the lo- 
cation (xi,yi). The smoothing kernel K is characterized 
by its smoothing bandwidth h, which we adopted as vari- 
able such as to include a fixed number n of stars within 
3h at each location {x,y). For convenience in comparing 
our analyses, we follow the prescription of Walker et al. 
(2006b) in adopting a bivariate Gaussian kernel; hence 

K DC exp|-| ^^~^'^'^^^~^*^' |. The number of nearest 

neighbors n was chosen sufficiently high to yield a statisti- 
cally significant estimate of each local dispersion measure- 
ment, but sufficiently low so that any real spatial informa- 
tion would not be averaged out. Hence we varied n from 
5 to 30, corresponding to at most about 25% of the entire 
sample size. 

The significance of any potential substructure in our 
data was then determined via detailed Monte Carlo mod- 
eling. For this purpose, 10'^ random data sets were gen- 
erated, where the {x,y) coordinates of our targets were 
preserved, but the observed velocities and error esti- 
mates were permuted with respect to their position. In 
this way, the sample mean and dispersion were retained, 
whereas any spatial information was dissolved. For each of 
these random data sets, the spatial dispersion distribution 
(eq. 3) and its maxima and minima were determined. The 
significance phot {Pcoid) of hot (cold) substructures at each 
spatial point was then defined by the fraction of random 
data sets with a maximum (minimum) dispersion, which 
is lower (higher) than the actual observed local dispersion 
a{x,y) (see Walker et al. 2006b). 

We do not find evidence of any localized cold substruc- 
ture with low radial velocity dispersions in our data (see 
Fig. 6). In both the GMOS + DEIMOS data sets of Leol 
and after inclusion of M98 targets, the significance of cold 
structures, Pcoid, does not exceed 50% for any choice of 
the number of neighboring points. The typical value of 
this significance lies at 30%. 

On the other hand, the presence of a physically real lo- 
cal maximum in the radial dispersion profile could indicate 
the presence of a locally hot structure at radii of -^3-4'. 
Although the probabilities for an occurrence of hot sub- 
structure increase towards the nominal tidal radius, where 
the data are only sparsely distributed, the formal signifi- 
cances Phot are generally well below the 80% level. Hence, 
the apparent localized kinematical structure is consistent 
with statistical fluctuations about a constant dispersion 
value at the Icr-level. In particular, there is no significant 
evidence of any localized structure at the radii in question 
- at loci between 3' and 4', the probability of an occur- 
rence of any dynamically hot substructure is less than 15% 
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for any reasonable choice of n. 

Since the submission of our paper, S06 and Bosler et 
al. (2006) have presented kinematic data sets for Leol. 
An analysis of the combined data sets (also including our 
data and those of M98) , which yields a sample of the or- 
der of 400 red giants velocities, will aid the investigation of 
the the full kinematical structure of Leo I and is currently 
underway. 

From our data, we conclude that the apparent rise in 
the velocity dispersion profile docs not cxclndc the pres- 
ence of any kinematical substructure, but it is probably 
not a localized feature. 

5.3. (N on-) Influence of binaries 

The presence of a significant population of binaries in 
any kinematical data set leads to an inflation of the ob- 
served velocity distribution, or, in other words, the true 
line of sight velocity dispersion of a stellar system is smaller 
than the observed dispersion, as soon as a non-zero binary 
fraction is considered. Moreover, the high-velocity tail of 
a binary distribution drives a distribution function away 
from Gaussian, so that an over-simple determination tends 
to produce larger errors in the derived velocity dispersion, 
which in turn may reduce the statistical significance of 
any real radial gradient in the dispersion profile (see error 
bars in Fig. 6). Fortunately, the effect of binary stars has 
been shown to be negligible in dSphs in the past. Harg- 
reaves, Gilmore & Annan (1996) could effectively show by 
using Monte Carlo simulations that the velocity disper- 
sion caused by pure binary orbits is small compared with 
the generally large observed dispersions in dSphs, provided 
that the ensemble of orbital parameters in the dSph stars 
is similar to the distribution in the solar neighbourhood. 
Likewise, Olszewski et al. (1996) concluded from their sim- 
ulations that none of the kinematical estimates in the UMi 
and Draco dSphs significantly changes under the influence 
of binary stars. This finding was underscored by repeat 
observations of red giants in the Draco dSph (Kleyna et al 
2002), and also Walker et al. (2006a) conclude from their 
repeat observation of red giant velocities in the Fornax 
dSph that the impact of binaries on the measured veloc- 
ity dispersion is negligible. These observed kinematics did 
not show any evidence that would support an overall bi- 
nary content larger than 40%. Moreover, the dynamically 
significant fraction in the Draco sample amounted to less 
than 5%. In the case of Leol, Gallart et al. (1999a) argued 
from their modeled SFH based on deep HST color mag- 
nitude diagrams (CMDs) that it is unlikely that the total 
binary fraction in this galaxy exceeds 60%, again with a 
far lower dynamically significant number. 

In order to explore the importance of binaries for the 
velocity dispersion profile of Leo I, we added an additional 
term to eq. 2, accounting for a binary star distribution 
Pb{v), which is then convolved with the observed veloc- 
ity distribution, in accordance with ccj. 1 of Kleyna et 
al. (2002). The binary probability distribution Pb{v) was 
adopted from Kleyna et al. (2002) and is based on the 
velocity measurements of Duquennoy & Mayor (1991) for 
a sample of solar-type primary stars in the solar neigh- 
bourhood. Kleyna et al. (2002) then obtained a realistic 
distribution by taking into account the giant branch binary 
evolution through circularization of the orbits. 



We note that the presence of a non-zero binary fraction 
does not alter the derived interloper fraction, which is still 
consistent with zero (less than 5% in the outermost bin). 
Fig. 7 shows velocity dispersion profiles obtained under the 
assumption of different binary fractions. The average er- 
rorbar on the velocity dispersion measurements increases 
by 10% for a binary fraction of 0.4 and by 25% for an 
implausibly large fb of 0.8. Likewise, the overall decrease 
of the dispersion profiles is, with less than 0.7kms~^ de- 
viation for 40% binaries, well below the 0.5 ct level. An 
increase of leads to a progressively more pronounced 
fall off of the underlying dispersion profiles. Nonetheless, 
the velocity dispersion profile assuming 40% binaries is 
consistent with that using no binaries at the 90% confi- 
dence level, as confirmed by a K-S test so that henceforth 
the effect of binaries on our observed line of sight velocity 
dispersion profiles is taken to be negligible. 

6. ISOTROPIC MASS ESTIMATES 

The simplest possible estimate of the mass profiles of 
dSphs uses single-component dynamical models, in which 
it is assumed that the mass distribution follows that of the 
visible component ("mass follows light", e.g., Richstone 
& Tremainc 1986). In this context, King models (King 
1966) were often used to describe both the surface density 
and the velocity dispersion profile. However, the parame- 
ters describing the visible content of dSphs, i.e., core ra- 
dius Tc and velocity dispersion profile cr, are in all cases 
studied in detail inconsistent with the mass-follows-light 
assumption inherent in this style of analysis (Kormendy 
& Freeman 2004). In the case here, where the approx- 
imately fiat observed dispersion profile of Leol suggests 
that the stars orbit in a (dark matter) halo which extends 
to radii larger than the nominal tidal radius of the observed 
light distribution, this historical approach is inappropri- 
ate. Where sufficient data exist, two parameter dynamical 
models for spherical stellar systems (Pryor & Kormendy 
1990; Kleyna et al. 2002), additionally accounting for ve- 
locity anisotropics and a dark matter component or non- 
parametric modelling schemes (Wang et al. 2005, Walker 
et al. 2006a) arc progressively being applied to reproduce 
observed profiles, in particular out to large radii. 

Where only limited data exist, sufficient to define a 
binned dispersion profile but insufficient to support a full 
dynamical analysis of the distribution function, an inter- 
mediate level of analysis is appropriate, and is applied 
here. We pursued a simple approach to obtain an esti- 
mate of the galaxy's mass and density profile by integrat- 
ing Jeans' equation (Binney fc Tremaine 1998, eqs. 4-54 
ff.) first under the assmnptions of an isotropic velocity dis- 
tribution and spherical symmetry, and below considering 
anisotropic distribution functions, using smooth functional 
fits to represent the light distribution and dispersion pro- 
file of Leo I. 

In the case of an isotropic velocity distribution, the 
Jeans equations give rise to the simple mass estimator 

««-^(^). 

where u and v'^ denote the deprojected, three dimensional 
light density distribution and radial velocity dispersions. 
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respectively. Both these quantities are obtained via di- 
rect deprojection of the observed surface brightness, I{r), 
and projected dispersion profile a{r) adopting a conve- 
nient functional form, which in this case is a Plummer 
model (Wilkinson ct al. 2006a, after Binncy & Trcmaine 
1987). In this case, the surface brightness and associated, 
deprojected 3D profiles read 

The surface brightness profile of Leo I, which we adopted 
from Irwin & Hatzidimitriou (1995), is well fit by a Plum- 
mer model (left panel of Fig. 8) with a scale radius a of 
4f6, which also nicely accounts for objects near the nom- 
inal tidal radius at 12.'6. Since the observed velocity dis- 
persion profile, on the other hand, appears to be consistent 
with being flat, we fit this component with a Plummer pro- 
file with a large scale radius of 74' to ensure a flat shape 
throughout the observed radius (right panels of Fig. 8). 
However, as Fig. 8 indicates, this argument is slightly sen- 
sitive to the value of the outer bin, which may give rise to 
a fall-off in the profile. Hence, we also chose a Plummer 
profile with a smaller radius (36'), which fits the observa- 
tions best. In either case, the central value of the velocity 
dispersion uo was determined as 9.8kms~^. These values 
hold only for the combined CMOS + DEIMOS sample. 
The inclusion of M98's data leads to a decreased central 
velocity dispersion and consequently yields a lower central 
density po- The resulting density profiles in Fig. 9 fall off 
faster in the case that Leo I has a falling dispersion profile. 
Either density profile reaches a r^^ slope at 3', but un- 
der the assumption of a flat velocity dispersion the density 
converges to the r~^-law (linearly rising mass), whereas it 
tends to be consistent with a j-"^ ''' behavior towards 
the nominal tidal radius for the falling dispersion profile. 

As discussed above our data suggest a rising feature 
in our velocity dispersion profile at approximately 3'. 
Although it has been shown above that this is not a 
statistically-significant localized structure, it may be phys- 
ically real, and so provide useful constraints on the dark- 
matter mass distribution in Leo I. Consequently, we con- 
sider this particular shape by fitting the dispersion profile 
with an additional asymmetric Gaussian component of the 
form (Jbump oc 10""'^ X exp(-10(''"^)/P) overlaid on the 
Plummer profiles. The free parameters a, /i and p are 
then determined in a least-squares sense. Using this best- 
fit representation in the dynamical calculations results in 
an unphysical behaviour in the resulting mass profile, i.e., 
a drop in the cumulative mass and the density profile at 
the location of the bump (left panels of Fig. 9). This is at- 
tributed to a too steep rise in the additional component as 
opposed to the smooth decline towards larger radii. Hence, 
the associated gradient, which enters the Jeans' equation, 
overwhelms the gradient of the declining portion in the 
dispersion and light profiles leading to a negative contri- 
bution in the derived mass profile. We note that varying 
the parameters of the Gaussian peak, retaining consistency 
with the GMOS + DEIMOS profile, does not remove this 
inconsistency. Conversely, a comparable fit to the GMOS 
-I- DEIMOS + M98 data does not exhibit this unphys- 
ical outcome, since for the respective best-fit profile the 



gradient in the innermost regions does not counteract the 
gradient of the more gradually declining profile at larger 
radii. We conclude that the failure of such a basic func- 
tional approach to simultaneously yield a physically expe- 
dient representation of all profiles considered suggests that 
we in fact do not see a real substructure. If, on the other 
hand, the dispersion profile could be modeled sufficiently 
well in either case and the radial feature was assumed to 
refiect a real structure, its particular shape would result 
in a close to uniform density profile at central radii. This 
would be indicative of a central core in Leo I. We will ex- 
plore the plausibility of cored mass profiles in the next 
Section. 

The upper and lower mass- and central density limits 
were finally derived by determining a profile that repro- 
duced the observed velocity dispersion profile within the 
measurement uncertainties on either side. An estimate 
of Leol's mass thus yields Mtot = (8 ± 2) x lO'^Af© en- 
closed within the King radius of 12f6, corresponding to 
0.9 kpc at the distance of Leo I. The exact values de- 
pend on the choice of the formal representation of the 
line-of-sight velocity dispersion profile and are detailed 
in Table 4. Similarly, the central density is found to be 
Po = (2.8 ± 0.8) X 10* Mq kpc-3, which is, in the light of 
the measurement uncertainties, consistent with the value 
of (3.4±0.9) X 10* Mq kpc"^ obtained by M98 from single- 
component King fitting. The available dynamical mass 
estimates for dSplis yield total masses out to the radial 
limits of the respective kinematic data in the range of 
3-8xlO'^M0 (e.g, Mateo 1998, Wilkinson et al. 2004; 
Kleyna et al. 2004; Chapman et al. 2005; Walker et 
al. 2006a; Wilkinson et al. 2006a). This places Leo I 
in the upper range of known LG dSph masses. Adopt- 
ing a total luminosity of L = (3.4 ± 1.1) x 10^ Lq (Irwin 
& Hatzidimitriou) we find an isotropic mass-to-light ra- 
tio of {M/L)tot = (24 ± 6) {M/L)q (see Table 4). This 
value is three times larger than the central value quoted 
by M98, who obtained their estimate based on analyses of 
the galaxy's core region and adopted a higher luminosity 
of L = 4.9 X 10^ Lq. If we adopt this latter value for L, we 
obtain a {M/L)tot of 17 {M/L)q. Moreover, S06 derive an 
estimate of 5.3 {M/L)q, again with an underlying higher 
total luminosity of 7.6 x 10^ Lq, and based on core fitting 
and their value of the central velocity dispersion. 

In Fig. 10 we plot the V-band mass to light ratios ver- 
sus V-band Imninosities of the known Galactic dSphs, in- 
cluding some of the recently discovered systems, like Ursa 
Major (Willman et al. 2005) and Bootes (Belokurov et al. 
2006a), plus two of the M31 companions with published 
masses. In particular, we adopted mass estimates from 
Mateo (1998) for Leo II and Sculptor, whereas more recent 
measurements were used for the other satellites (Cote et 
al. 1999 for And II; Wilkinson et al. 2004 for UMi and 
Draco; Kleyna et al. 2006 for Ursa Major; Chapman et al. 
2005 for And IX; Wang et al. 2005 for Fornax; Wilkinson 
ct al. 2006a for Carina and Sextans, and Munoz et al. 
2006b for Bootes). 

The idea that dSphs may be dominated by dark matter 
out to large radii raises the intriguing question whether all 
these dwarf galaxies could be enclosed in comparable dark 
matter halos of similar total mass, as first pointed out by 
Mateo et al. (1993). This feature would then argue in fa- 
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vor of an intimate relation between the dark matter halos 
of the dSph systems. Moreover, cosmological simulations 
suggest that it is possible that there is a minimum mass 
for a dark matter halo that contains stars (e.g., Read et 
al. 2GG6a and references therein). It is also an interesting 
question to ask whether all the dSphs have similar mass 
haloes, as this might tell us something about the proper- 
ties of the dark matter. Under such an assumption, the 
halos can be represented by the the simple relation 

{M/L)tot = [M/L), + Mdm/L, (6) 
where {M/L)^, is the intrinsic mass-to-light ratio of the 
stellar component in the galaxy, Mdm denotes the mass 
of the dark matter halo and L is the integrated luminosity 
in the V-band (M98; Wilkinson et al. 2006a). It is realistic 
to assume that the (M/L)-ratio of the the luminous com- 
ponent in dSphs is of the order of the mean value observed 
in low-concentration globular clusters, which are charac- 
terised by low values implying no dark matter (Pryor et 
al. 1989; Moore 1996; Dubath & Grillmair 1997; M98; 
Baumgardt et al. 2005). We do not correct the observed 
M/L-ratios for the (small) effects of stellar evolution, ac- 
counting for each individual dSph's SFHs. Thus we set 
(M/L)* = 1.5(M/L)q. As the best fit dashed line in 
Fig. 10 implies, the observed dSphs scatter around the ex- 
pected relation for a population of luminous objects that 
are embedded in a dark halo with a common mass scale of 
the order of 3 x 10^ M©. Considering the generally large 
scatter in such a (M/L)-My-diagram, Leo I fits well into 
the global picture and apparently is governed by the same 
uniform dark halo properties. It is hence worth noticing 
that Fig. 10 demonstrates that the narrow range of dSph 
velocity dispersions of 6-10 km s~^, combined with their 
rather similar physical length scales, can generally be in- 
terpreted in terms of a common mass scale so that this 
kinematical piece of information can reliably be used as a 
proxy for an overall underlying mass distribution (see also 
Wilkinson et al. 2006a). 

7. VELOCITY ANISOTROPY 

The mass estimates derived above are based on the as- 
sumption of an isotropic velocity distribution, i.e., the ve- 
locity anisotropy parameter (3 = 1 — {vg)/(vf) was as- 
sumed to be zero. However, most of the kinematic stud- 
ies of dSphs have revealed that a non-negligible amount 
of anisotropy is required to account for the shapes of the 
galaxies, and the observed dispersion profiles (Lokas 2001, 
2002; Kleynaetal. 2002; Wilkinson et al. 2004). Since the 
neglect of a non-zero anisotropy has the effect of under- or 
overestimating the dark halo mass, we also shall consider 
this problem by solving the Jeans equation 

-u ^ (7) 



(1 r r 



dr 



for a varying (3. Here, v and denote again the 3- 
dimensional light and dispersion profiles, and $ is the grav- 
itational potential associated with the underlying mass 
distribution M(r). The actual observable quantity, the 
line-of-sight velocity dispersion a, is then obtained by di- 
rect integration along the line of sight, and under the con- 
straints pv^ ^ for r — > oo and /3 = const. This gives 
rises to the one-dimensional expression 



-20+1 



(8) 



(Binney & Mamon 1982; Lokas & Mamon 2003; Wilkin- 
son et al. 2004). By adopting analytical prescriptions for 
the involved functions, one can then proceed to derive a 
numerically. 

In concordance with the previous section, we chose to 
describe Leol's surface brightness profile I{R) and the as- 
sociated 3D-deprojection with a Plummer-profile using the 
best-fit parameters to the observations. We then use two 
different representations of the mass profile M(r), which 
enters this formalism, spanning the range of plausible mass 
distributions. 

7.1. NFW halo mass distribution 

Based on high-resolution cosmological simulations, 
Navarro, Frenk & White (1995, hereinafter NFW) demon- 
strated that the density profiles of dark matter halos are 
well fit by a simple function with a single free parameter, 
the characteristic density. This solution has been applied 
to a wide range of halo masses, ranging from galaxy clus- 
ters to small scales such as the cosmological dark matter 
halos associated with the dSphs. Upon transformation of 
the variables, the density profile corresponds to a mass 
profile of the form 

where s = r/r-^ denotes the radial distance in units of the 
virial radius and My is the mass enclosed within the virial 
radius. The latter is generally identified with the total 
mass of the halo. Finally, the concentration parameter c 
is used to describe the shape of the profile and defines the 
amplitude function 

.g(c) = [ln(l + c) - c/(l + c)]-\ (10) 

For a detailed discussion of these parametrizations and the 
model, we refer the reader to NFW or Lokas & Mamon 
(2001). Since the NFW-profile has the unphysical disad- 
vantage of diverging at large radii, it is convenient to con- 
sider only radii within a certain cut-off radius. We follow 
long-standing practice and assume this cut-off to coincide 
with r^; beyond this point, the overall density distribu- 
tion becomes imreliable in most natural cases (e.g., Lokas 
& Mamon 2001). 

The concentration c has been shown to scale with mass. 
Here we adopt an extrapolation of the formulae derived in 
the A''-body simulations for ACDM cosmology of Jing & 
Suto (2000) to the small masses of dwarf galaxies. Hence, 

where h denotes the Hubble constant in units of 
lOOkms"^ Mpc~^. In concordance with current cosmo- 
logical results, we will use h=0.7 for the remainder of this 
work. Likewise, the virial radius is related to the virial 
mass via 
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1/3 



(12) 



^1012Mq, 

This leaves My and the anisotropy /3 (assumed to be con- 
stant) as the only free parameters to be determined in a 
fit of the NFW model profile to the observations. 
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As Fig. 11 (top panel) implies, none of the values for the 
anisotropy which wc have considered yields an adequate fit 
of the overall dispersion profile. It rather appears that a 
progressively radial anisotropy is needed to explain the ob- 
servations at larger radii, whereas the inner parts of the ve- 
locity dispersion are best represented by an isotropic veloc- 
ity tensor. It is worth noticing that the innermost bin, with 
a smaller value for the dispersion in the sample after inclu- 
sion of the M98 data points, may suggest the presence of 
some tangential anisotropy in the innermost regions. How- 
ever, based on the observations, /3 is unlikely to fall below 
—0.5. For the four different (3 curves that we tested (i.e., 
—0.5, 0, 0.5 and 1), wc obtained reduced values of (0.2, 
0.16, 0.22, 0.78). The final best fit in terms of a minimized 

was achieved for (3 = 0.05 so that the overall profiles are 
in agreement with an isotropic velocity distribution. The 
resultant best-fit total virial mass is My = 1.0 xlO^M©, 
with a corresponding virial radius of r„ = 21 kpc and a 
concentration c — 19.4. Furthermore, the mass within the 
nominal observed tidal radius amounts to 7.7x10"^ M©. It 
is, unsurprisingly given the near- isotropic distribution de- 
duced, in good agreement with the value obtained under 
the assumption of a purely isotropic velocity distribution. 

7.2. Cored halo profile 

The density profile that we obtained under the assump- 
tion of velocity isotropy appears to be indicative of a flat, 
close to uniform density core at small radii (see bottom 
panels of Fig. 9). In conjunction with the observational in- 
dications of cored density profiles in low-luminosity galax- 
ies and in particular in some of the dSphs (Lokas 2002; 
Kleyna et al. 2003; Strigari et al. 2006; Wilkinson et al. 
2006a) this prompted us to also consider this particular 
density distribution in our calculations of Leol's disper- 
sion profile (eq. 7). We parameterize the cored density 
profile in terms of a generalized Hernquist profile (Hern- 
quist 1990; Zhao 1996; Read & Gilmore 2005) as 



p(r) = C 



-7 r 



1 



(7-3)/a 



(13) 



where C is a normalization constant, Vs is a scale radius, 7 
is the log-slope of the density profile within and a deter- 
mines the smoothness of the transition towards the density 
profile at larger radii. The case of a = 1 has often been 
discussed and fit to dispersion profiles in the literature, 
since it allows one to envisage a broad class of dark mat- 
ter halo representations (Jing & Suto 2000; Lokas 2002). A 
slope of 7 > then corresponds to the cuspy halos, where 
NFW-like profiles as discussed above are represented by 
7=1. The opposing class of models with a cored halo is 
realized by 7 = 0. 

It turns out that a least squares fit under the constraint 
a = 1, yields an unsatisfactory representation of our ob- 
servations. Thus we chose to use the profile according to 
eq. 13 in our fits, leaving C, rg, a and 7 as free param- 
eters. The resulting parameters of = 0.3 kpc, a = 1.6 
and 7 = improved the reduced from 55 (3-parameter 
fit) to 13 (4 parameters). It is worth noticing that indeed 
the best fit was obtained for a pure core profile (7 = 0). 
Fig. 11 (bottom panel) displays the resulting radial ve- 
locity dispersion profiles for varying degrees of (constant) 
anisotropy. 



As for the case of the cuspy dark halo treated above, we 
find that our observed velocity dispersion is broadly consis- 
tent with an isotropic velocity distribution,where the best 
fit yielded a /3 of 0.05. Still, the observations allow for a 
wider range in the anisotropy, from —0.5 to 0.5. In this 
context, a slight amount of tangential anisotropy can ac- 
count for the shape of the profile toward the inner regions, 
whereas significant radial anisotropy is ruled out at all 
radii. The values for the associated reduced x^ statistics 
for the cored halo are (0.39, 0.25, 0.40, 1.94) for respective 
j3 values of (-0.5, 0, 0.5, 1). 

Finally, we also adopted a radially varying anisotropy in 
our calculations, where we parameterized /? according to 
the traditional prescriptions by Osipkov (1979) and Mer- 
rit (1985) as (3oM{r) = / (r^ -I- r^), where Va denotes the 
anisotropy radius, at which the transition from isotropic 
to radial orbits occurs (see also Lokas 2002). Under the 
assumption of a cored density profile and this particular 
P{r), we can fit our observed line-of-sight velocity disper- 
sion profile fairly well (see the thick solid line in Fig. 11, 
bottom panel). The best-fit (at a x^ of 0.02) is then ob- 
tained for an anisotropy radius of Tq = 0.13 kpc. In partic- 
ular, the resulting curve is able to account for the observed 
rising shape of the profile discussed in Sect. 5.2, albeit by 
requiring a considerable amount of radial anisotropy at 
relatively small radii. 

Another possible explanation of the dispersion profile 
of Leo I might be the presence of two stellar components 
with different length-scales and associated different veloc- 
ity dispersions. McConnachie et al. (2006) have recently 
presented dispersion profiles derived for such a system. 
Their profiles bear a striking resemblance to that of Leo I 
and might account for the initial rise in the projected dis- 
persion. 

In concluding we note that the enclosed mass within the 
tidal radius amoimts to 8.5 x 10^ Mq under the assumption 
of a cored halo in Leo I. This value is larger than under the 
assumption of a NFW-like mass distribution, but still in 
agreement with the numbers derived from the Jeans equa- 
tions above. 

8. THE METALLICITY OF LEO I 

The observed wavelength region in the near-infrared and 
our spectral resolution, coupled with the S /N of our obser- 
vations also enabled us to get an estimate of the red giants' 
metallicites from the equivalent widths (EWs) of the Ca 
triplet (CaT) . This procedure succeeded for 58 of the radial 
velocity member stars targeted with GMOS, all of which 
lie well within the nominal tidal radius. For the remainder 
of our targets, the S/N did not allow us to reliably mea- 
sure EWs, where a lower limit for these measurements was 
S/N~ 10. 

8.1. Calibration of the metallicity scale 

The prominent Ca II feature may be calibrated onto 
global metallicities since the linestrength is a linear func- 
tion of metallicity for red giants. By defining a reduced 
width W , which accounts for the luminosity difference be- 
tween the targeted star and the horizontal branch (HB) of 
the system, one can effectively remove any dependence on 
stellar gravity, effective temperature, and distance. Thus, 
the CaT has become a well defined calibrator for assess- 
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ing the nictallicity for old, globular cluster-like populations 
(Armandroff & DaCosta 1990, Rutledge et al. 1997a,b). 
Consequently, the derived line widths are generally cali- 
brated onto reference scales for globular clusters of known 
metallicities (Zinn & West 1984; Carrctta & Gratton 1997; 
Kraft & Ivans 2003). Leo I has a prominent old population 
(e.g.. Held et al. 2000), but its dominant populations are 
of intermediate age (e.g., Gallart et al. 1999a, b). Nonethe- 
less, we can still apply the CaT technique, since Cole et 
al. (2004) have extended its calibration to younger ages. 
These authors demonstrated that the CaT technique can 
be used within an age range from 2.5 to 13 Gyr and for 
a metallicity range spanning —2 < [Fe/H] < —0.2. Since 
our Leo I observations aimed primarily at measuring accu- 
rate radial velocities, we did not target any globular clus- 
ters, which could have been used as calibration standards 
for the CaT technique. Hence, we had to rely on stan- 
dard calibrations devised in the literature. Following long 
standing practice, we employed the definition of Rutledge 
et al. (1997a, hereafter R97a) for the linestrength of the 
CaT as the weighted sum of the EWs, giving lower weight 
to the weaker lines: 

SVF = 0.5W8498 + W8542 + 0.6W8662. (M) 

For a few of the spectra the Ca line at 8498A was too weak 
to be accurately measured or was hampered by too low a 
S/N over the respective bandpass. In this case we used a 
relation derived from the high S /N spectral measurements 
of a sample of Galactic globular clusters from Koch et al. 
(2006a): 

SVF = 1.13(W8542 + 0.6W8662) + 0.04. (15) 

This relation was applied to the low-quality spectra where 
only the two strongest lines had S/N>10. R97a defined 
the reduced width as 

W' = i:W + 0.64 (±0.02) {V - Vhb) , 

where TW denotes the CaT linestrength and \ hb is the 
magnitude of the horizontal branch of Leo I. For Leo I, the 
HB luminosity is (22.60 ± 0.12) mag according Held et al. 
(2001), who measured it based on this galaxy's RRLyrae 
population. 

The final calibration of the reduced width in terms of 
stellar metallicity is then obtained via the linear relation 

[FelH]cQ = -2.66 (±0.08) -|- 0.42 (±0.02)VF' 

(Rutledge et al. 1997b), where we chose to tic our obser- 
vations to the reference scale of Carretta & Gratton 1997, 
hereafter CG97) unless stated otherwise. 

The EWs of the calcium lines were determined in anal- 
ogy to the methods described in Koch et al. (2006a), i.e., 
by fitting a Gaussian plus a Lorentzian profile across the 
bandpasses defined in Armandroff & Zinn (1990) and sum- 
ming up the flux in the theoretical profile. The median 
formal measurement uncertainty achieved in this way, in- 
corporating EW measurement and calibration errors, is 
0.12 dex. However, we note that the intrinsic accuracy of 
the CaT based metallicity measurements may be signifi- 
cantly lower. For instance, variations in the HB level of a 
composite stellar system with both age and metallicity can 
give rise to a systematic error of the order of 0.05-0.1 dex 
(Cole et al. 2004). Moreover, the general incompatibilitiy 

* We note that this metal poor regime was not sampled by the GCs 
extrapolation of the existing calibrations. 



between the [Ca/Fe] ratio in dSphs and the Galactic glob- 
ular clusters adopted as calibrators is another source of 
uncertainty. High resolution abundance studies of red gi- 
ants in dSphs have shown that their a-element abundance 
ratios lie significantly below those in Galactic halo stars 
at comparable metallicities (e.g., Venn et al. 2004), and 
also the two red giants observed in Leo I by Shetrone et 
al. (2003) follow this trend of a-deficiency. Moreover, the 
low-resolution studies of Hosier et al. (2006) underscore 
these low abundance ratios as c;ompared to the halo. The 
systematic uncertainties resulting from such a priori un- 
known variations in [Ca/Fe] for our targets can reach up 
to 0.2 dex (e.g., Koch et al. 2006a). It is worth noticing 
that such effects can be overcome in future studies by cal- 
ibrating the reduced width directly onto [Ca/H] (Hosier et 
al. 2006), thus circumventing any dependence of the cali- 
brations of the galaxy's SFH, in contrast to the presently 
employed method. 

8.2. Metallicity distribution 

Fig. 12 shows the resulting histogram of our metallic- 
ity estimates. The MDF is clearly peaked at a median 
metallicity of —1.31 ± 0.02 dex on the scale of CG97, with 
a measurement error of 0.12dex (—1.61 on the scale of 
Zinn & West 1984; hereafter ZW84). This value is in ex- 
cellent agreement with the study of Hosier et al. (2006), 
who derived a CaT metallicitity of — 1.34±0.02dex (mean 
measurement error of 0.10) for a sample of 101 red giants 
(note, however that the lack of any common targets be- 
tween the studies of Hosier et al. [2006] and ours does not 
allow us to compare measurements of individual stars). 
Our results also agree within the uncertainties with the 
low-resolution studies of Suntzeff et al. (1986), who ob- 
tained — 1.9 ± 0.2 dex (ZW84). Moreover they are consis- 
tent with the estimates from photometric analyses, which 
broadly agree on a mean metallicity around —1.6 ±0.4 dex 
(ZW84 scale; Demers et al. 1994; BcUazzini ct al. 2004) 

The formal 1 cr-width of the MDF is 0.25 dex. Taking 
into account a broadening due to the observational un- 
certainties, the intrinsic metallicity dispersion amounts to 
0.22 dex. With a full range in [Fe/H]cG97 of approximately 
one dex between —1.8 and —0.8 dex (note the most metal 
poor star, which we discuss below) we observe a smaller 
spread than is seen in a number of other dSphs (Shetrone, 
Cote, & Sargent 2001; Pont et al. 2004; Tolstoy et al. 2004; 
Koch et al. 2006a), where the MDF may span more than 
2 dex and exhibits a well populated metal poor tail. Leo I 
has been shown to exhibit a rather narrow RGB. By means 
of detailed modeling of Leo I's CMD and in particular of 
this narrow feature, Gallart et al. (1999a) suggested that 
this might be indicative of the absence of any significant 
chemical enrichment in this dSph. Clearly, the spectro- 
scopic data show that this is not correct and that Leo I 
underwent considerable enrichment of at least 1 dex in 
total range. 

There is one object in the sample, #774, for which we 
derived a CaT based metallicity of —2.61 dex. There is 
no evidence in the spectrum that the linestrength of this 
target has been underestimated, since its S/N of ^45 al- 
lowed an accurate determination of the triplet's EWs (see 

in R97b so that the determination of such metallicities relies on an 
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Fig. ISjfl. Although weak neutral metal lines around the 
CaT region do not show up clearly at our nominal spectral 
resolution and are thus mostly blended, there are still some 
additional absorption features discernible from the noise 
(see Fig. 13). Among these are for instance the distinct 
Fe I blend at ~8514A and a visible Fe I-Hne at -8689A. 
Though clearly present in our stars with highest metallici- 
ties, these features are hardly detectable in the continuum 
of the metal poor candidate. By assuming appropriate line 
bandpasses around these two iron features, we employed 
the same fitting technique as for the CaT described in 
Section 8.1 in order to derive a rough qualitative estimate 
of these lines' EWs in each of the stars, where possible. 
From this point of view, #774 is well consistent with be- 
ing a rather metal poor object, in which these two EWs 
are practically governed by the continuum noise. Further- 
more, its location well within the galactic boundary (2.5 
core radii, 0.67 tidal radii, resp.), its consistency with the 
systemic velocity (99.9% confidence level) and the low in- 
terloper fraction around Leo I's radial velocity make it ap- 
pear unlikely that this star is a contaminating foreground 
dwarf. 

8.3. Radial gradients and implications for evolutionary 

models 

In those dSphs of the Local Group, in which both old 
and intermediate-age stellar populations have been de- 
tected, the younger components often are more strongly 
centrally concentrated. In a number of dSphs with pre- 
dominantly old populations, population gradients have 
been found in the sense that red HB stars are more cen- 
trally concentrated than blue HB stars (e.g., Da Costa et 
al. 1996; Hurley-Keller et al. 1999; Harbeck et al. 2001; 
Tolstoy et al. 2004), a trend that is mirrored also by the 
presumably more metal-rich red giants. However, as also 
pointed out by Harbeck et al. (2001), this trend is not seen 
in all dSphs. 

Since Leo I hosts dominant intermediate-age popula- 
tions as well as old populations, one may expect that the 
more metal-rich and presumably younger stars should then 
also be more concentrated with respect to a spatially ex- 
tended metal-poor population. Our data suggest, how- 
ever, that the metal-rich and metal-poor components in 
our MDF are drawn from the same spatial distribution. 
This is reflected by the comparison of the cumulative radial 
distributions of the more metal-poor ([FeH] < 1.3 dex) ver- 
sus the metal-rich ([FeH] > 1.3 dex) populations in Fig. 14. 
Via a K-S test we could underscore this lack of a radial 
metallicity gradient at the 93% confidence level. This find- 
ing is in concordance with the lack of a considerable pop- 
ulation gradient in Leo I's CMD. As was reported by Held 
et al. (2000), both the old HB population and its numerous 
intermediate-age counterpart exhibit essentially the same 
spatial distribution (see also the right bottom panel of 
Fig. 14). In this respect, Leol differs from other dSphs 
with prominent intermediate-age populations such as Ca- 
rina (e.g., Harbeck et al. 2001; Koch et al. 2006a), Fornax 
(Grebel 1997; Stetson et al. 1998) and Leo II (Bellazzini 
et al. 2005). Moreover, we find that the kinematics of the 
metal poor and the more metal rich samples do not differ 
significantly, neither. This is also in contrast to what is 
found for different stellar components in many other dSphs 



(Tolstoy et al. 2004; Ibata et al. 2006): While the sample 
with [Fe/H]< —1.3 exhibits a dispersion of 9. l±2.8kms~^, 
for those stars with [Fc/H]> -1.3, a = 10.7±3.1 kms^^, 
which is an insignificant difference in the light of the un- 
certainties. 

This behavior may indicate that Leo I's star formation 
occurred from a well mixed reservoir of gas that was little 
affected by local accretion of material or spatially local- 
ized outflows. In this vein, it was suggested by Held et 
al. (2000) that the size of this galaxy has not significantly 
changed and that it is unlikely to have undergone severe 
structural changes in the course of its evolution, i.e., since 
the first onset of star formation. In such a model one 
must assume that external effects such as tides, accretion 
or minor mergers may not have played a significant role in 
Leo I's history, neither kinematically nor chemically. On 
the other hand, S06 claim that the preferential periods of 
SF in Leo I may be linked to perigalactic passages, which 
argues in favor of tidal interactions. Still, without extant 
detailed knowledge of this dSph's orbit (via proper mo- 
tions), this has to await further investigation. 

8.4. Comparison with simple chemical evolution models 

Present-day chemical evolutionary models agree on the 
fact that dSphs are mainly controlled by low SF efficiencies 
compared to the solar neighborhood, and by strong galac- 
tic winds of the order ten times the SF rate (Lanfranchi 
& Matteucci 2004). For want of a detailed model predic- 
tion of Leo I's MDF from its SFH we follow long standing 
practice in applying basic models of chemical evolution to 
describe the observations (e.g., Pagel 1997, ch. 8), without 
employing any deeper nucleosynthetic or physical prescrip- 
tions. Such more detailed modeling is left for future work. 

In that spirit, Fig. 12 shows the predicted MDFs for 
two simple models of chemical evolution. To start with, 
we obtained a best-fit representation of a modified sim- 
ple closed-box model allowing for outflows, which is exclu- 
sively parameterized by the effective yield Poff- Basically, 
this value is proportional to the true nucleosynthetic yield 
and takes into account effects of the loss of metals. Conse- 
quently, the mean metallicity is reduced to below the true 
yield by the same factor, which relates the outflows and 
the actual SF rate. A low yield of 0.06 is well able to 
reproduce the metal poor mean location of the peak of the 
MDF, when compared to the solar neighborhood, which 
peaks at around —0.2 dex (Nordstrom et al. 2004 and ref- 
erences therein). However, this prediction tends to over- 
estimate the number of stars below —1.8 dex (note that 
the most metal poor object in our sample was excluded 
from the fit) , leading to the well-known G-dwarf problem. 
We note that a comparison of these models, which are 
based on predictions for long-lived stars, to our observa- 
tions of K-giants with predominantly negligible lifetimes 
is still valid, since the simple model of chemical evolution 
does qualitatively predict an excess of metal poor K-giants 
under the assumptions of a standard IMF and loss of met- 
als from the galaxy (Koch et al. 2006b). Moreover, it has 
been shown that there is no noticeable difference between 
the K-giant MDF in the solar neighbourhood (Mc William 
1990) and the local G-dwarf MDF from, e.g., Nordstrom 
et al. (2004). Hence we carry on with the long-standing 
practice of applying the G-dwarf predictions to K-giant 
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MDFs as tracers of low-mass stars. 

After allowing for an admixture of a non-zero initial 
metallicity (prompt initial enrichment, PIE), the overpre- 
diction toward the metal poor tail is reduced and it turns 
out that a prc-cnrichmcnt of Zq = 0.013 Zq yields a re- 
markably good representation of the overall shape of the 
observed MDF. We refer the interested reader to Koch et 
al. (2006a, 2006b) for a discussion of these standard mod- 
els. 

9. SUMMARY AND CONCLUSIONS 

We have obtained radial velocity measurements for 120 
red giants in the Leo I dSph. These data not only increase 
existing published samples by more than a factor of three, 
but also extend in radial distance beyond previous studies 
of the central regions of the galaxy (M98). 

Although the weak velocity gradient that our data ex- 
hibit coincides with an apparent rotation about the ob- 
served major axis of Leo I, the associated amplitude has 
been shown to be insignificant. Moreover, we do not detect 
any trend of mean radial velocity with major axis distance, 
and in particular we do not find any evidence of an asym- 
metric velocity distribution towards the outer regions, in 
the sense that the eastern (western) halves of the galaxy 
would preferentially exhibit an excess of low (high) velocity 
outliers outside the formal tidal boundary. Supplemented 
by Leol's high radial velocity, which we confirmed, and 
its large distance from the Milky Way, Leo I most likely 
represents an isolated system, which is currently not af- 
fected by Galactic tides, at least from a kinematical point 
of view. This is also consistent with our detection of nine 
radial velocity members outside the nominal 'tidal' radius, 
and implies that the difficulty of the King-model func- 
tional fit to surface brightness data at some (outer) ra- 
dius is not necessarily a dynamical evidence of a physical 
effect. 'Tidal' radii for dSphs, as typically derived from 
photometry, are labels, and are not proven to be related 
in any way to gravitational tides. The results of S06 indi- 
cate that Leol's surface brightness profile exhibits multi- 
ple breaks at large radii, where it deviates from common 
King models, while further recent photometric analyses do 
not provide any such clear evidence of photometric com- 
plexity above the noise level (Smolcic et al. 2006). Still 
the detection of dSphs' stellar populations far out beyond 
the limiting radius, as in the case of Carina (Munoz et al. 
2006a), suggests that some dSphs may in fact be perturbed 
systems, at least in their outermost regions. Clearly this 
has implications for estimates of the dark matter content 
of dSphs at the largest radii, as tidal effects may inflate 
the projected dispersion in these regions. 

Our derived projected radial velocity dispersion profile is 
flat out to the nominal tidal radius, a result which is turn- 
ing out to be a very common feature among the dSphs of 
the MW system, being established in numerous studies of 
large radial velocity samples out to and even beyond their 
nominal tidal radii (Mateo 1997; Kleyna et al. 2001, 2003, 
2004; Lokas 2001, 2002; Wilkinson et al. 2004; 2006a,b; 
Mufioz et al. 2005; Walker et al. 2006a; Westfall et al. 
2006). Coupled with the overall high (M/L) ratios, and 
lack of unambiguous evidence for non-equilibrium effects, 
the prime explanation of flat dispersion profiles is that 
dSphs are dark matter supported at all radii, although the 



N-body simulations of S06 indicate that also tides could 
reproduce such features in Leo I. 

Nevertheless, at its large present Galactocentric dis- 
tance, tidal influences would be plausible only if Leol's 
orbit were close to radial. However, our current data set 
does not allow us to constrain its proper motions and any 
information about its orbit will have to await future accu- 
rate proper motion measurements. 

It turns out that the observed dispersion profiles are 
largc^ly c;onsistent with an isotropic velocity distribu- 
tion, but non-negligible amounts of radial (tangential) 
anisotropy may be needed to account for the shape of the 
dispersion profile at the largest (smallest) distances. 

An apparent rise in the velocity dispersion profile at 
nearly the core radius was demonstrated not to be a lo- 
calized dynamical substructure in the galaxy. By using 
simple prescriptions to describe the shape of the observed 
velocity dispersion profile we obtain an estimate of the 
enclosed mass within the nominal 'tidal' radius of 900 pc 
of ^8x10"^ Mq under the assumption of velocity isotropy. 
Such a mass and the resultant high mass-to-light ratio 
of ~24(M/L)0 are higher than previously derived values 
(M98; S06). From their (M/L) of -6 13 (M/L)0, M98 
concluded that Leo I contains a significant dark compo- 
nent. These estimates are in reassuringly good agreement 
with the values derived from light profiles obtained from 
high-quality photometric data under the assumption of an 
extended constant dark matter halo density (Smolcic et 
al. 2006). Hence, our results strengthen this view and are 
in concordance with an extended dark matter distribution 
with a halo mass of a typical scale (~3xlO'^M0), which is 
shared by the other dSphs of the LG. 

By adopting a NFW density profile for Leo I, we esti- 
mate its total virial mass within its dark matter halo to 
be ~1 xlO^M0. The mass within the tidal radius is esti- 
mated to be ~7.7 xlO^Af0. This docs not change much 
under the assumption of a constant density core in Leo I, 
where the enclosed mass and best-fit anisotropy are con- 
sistent with the values from a cusped profile. Although 
the density profile obtained using Jeans equations, under 
the assumption of velocity isotropy, indicates a flat central 
core, this scenario is not significantly preferred against the 
NFW case if one exclusively argues on the base of the ve- 
locity dispersion profile. In general, the shape of the ob- 
servations seems to indicate that dark matter in fact domi- 
nates this dSph at large and small distances. We are aware 
that analysis only of a dispersion profile a{r) is not suffi- 
cient to derive a unique kinematical distribution function 
of a stellar system. Unless more numerous data sets and 
in particular the higher moments of the velocity distribu- 
tion are considered, there remains a degenerac;y between 
anisotropy and density, leading to identical line-of-sight 
dispersion profiles (e.g., Merrifield & Kent 1990; Lokas & 
Mamon 2003). While in general such an analysis requires 
more data than is currently available for Leo I, this will 
be pursued in future papers dealing with more detailed 
dynamical modeling (Wilkinson et al., in preparation). 

We derived metallicites for 58 of the member stars from 
the well established CaT calibration. In compliance with 
previous estimates that were drawn from both photomet- 
ric and low-number spectroscopic analyses we find a mildly 
metal poor mean of — l.Sldex (CG97 scale) and a full 



Stellar Kinematics and Metallicities in Leo I 



15 



spread in metallicity covering 1 dex. This is noticeably 
smaller than the ranges found in dSphs with extended or 
complex SFHs (Grebel 1997; Pont et al. 2004; Tolstoy et 
al. 2004; Koch et al. 2006a). This spread exceeds the pre- 
dictions from CMD modeling if one accounts for the nar- 
row RGB and Leo I's existing intermediate-age population. 
Hence, there may be also a large range in ages present in 
Leo I, as hinted at by the presence of also an old population 
(Held et al. 2000, 2001) and the preliminary age determi- 
nations by Hosier et al. (2004). Practically all dSphs that 
contain significant intermediate-age populations do show a 
central concentration of their young components (Stetson 
et al. 1998; Harbeck et al. 2001; Bellazzini et al. 2005). 
Intriguingly, Leo I's prominently younger and presumably 
metal rich population does not exhibit this distinct pop- 
ulation gradient. Neither its stellar populations (Held et 
al. 2000) nor the metallicities we derived here show any 
evidence of radial gradients. In this respect, Leo I is rather 
unusual among these systems. This property may again 
point to a negligible role of external influences, such as 
localized accretion of material, in Leo I's evolution, allow- 
ing for good mixing of its gas. Since tides would affect 
this mixing, it is therefore possible that its current large 
distance from the MW may be representative of its whole 
evolution and that it never approached the Galaxy very 
closely. This view is consistent with the finding that the 
fraction of intermediate-age populations is closely related 
to the dSphs' Galactocentric distance (e.g., Grebel 2003 
and references therein). This behavior can be taken as 
indicative of an increased influence of (e.g., ram pressure) 
stripping close to the MW (van den Bergh 1994; Grebel 
1997; Grebel et al. 2003) and will not have affected Leo I to 
a major extent. We note, however, that SOB argue that the 
presence of the various stellar populations in Leo I may in 
fact be linked to much closer perigalactic passages, which 
facilitated the onset of SF, a scenario that would then in- 
corporate a non-negligible role of tides in Leo's evolution. 

We identified one comparatively metal poor star in our 
sample, which is unlikely to be a foreground contaminant. 
With its likely [Fe/H] of — 2.6 (CG97 scale) it lies 0.8 dex 
below the metal poor tail of our observed MDF. It has 
been suggested that the explosion of only one massive pre- 
galactic Population III star produces a sufficient amount of 
energy to expel the entire gas from the first (high-z) cos- 
mological minihalos (Bromm et al. 2003) . If the expelled 
gas falls back afterwards, the subsequently forming gener- 
ation of stars, which we would observe as the most metal 
poor and the oldest stars in dSphs, should bear the chem- 
ical signature of these Population HI stars (e.g., Beasley 
et al. 2003; Kawata et al. 2006). This makes metal-poor 
stars in dSph galaxies the prime targets for searches of 
the signatures of the first stars. However, no star with a 



metallicity below ~ —3 dex has been detected in a dSph 
so far (Fulbright, Rich & Castro 2004; Sadakane et al. 
2004) . The question of whether the lack of any such ex- 
tremely metal poor objects in the current data is merely 
due to the low-number statistics of present-day observa- 
tions, hence reflecting an incomplete sampling, or whether 
such objects generally do not exist in dSphs remains un- 
solved. Thus the presence of such a comparatively low 
metallicity object in our data provides an invaluable test- 
ing ground to study the early phases of general galactic 
evolution. 
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2004 12 12 


CMOS 


3x1800 








2004 12 15 


CMOS 


3x1800 


Field 3 


10 08 39 


+12 19 43 


2005 01 06 


CMOS 


2x1800 


#1 


10 09 14 


+12 18 00 


2005 02 06 


DEIMOS 


7200 


#2 


10 07 40 


+12 17 00 


2005 02 07 


DEIMOS 


9000 


#3 


10 08 05 


+12 18 00 


2005 02 08 


DEIMOS 


5400 



Table 2 

Characteristics of the target stars. Note that the velocities from both the CMOS and DEIMOS samples have 

BEEN SHIFTED TO A COMMON MEAN. 



ID 


a 


S 


I 


V-I 


r 




Svhc 


R 


S/N 




f> [i"e/H\ca 




(J2000) 






['] 


[km s 












103 


10 08 28 


12 18 51 


17.77 


1.43 


0.68 


288.06 


3.99 


59.78 


50 


-0.99 


0.13 


132 


10 08 48 


12 19 44 


18.03 


1.48 


5.15 


282.07 


6.28 


36.12 


28 


-1.45 


0.11 


137 


10 08 27 


12 20 31 


17.93 


1.56 


2.80 


292.81 


5.33 


43.99 


40 


-1.38 


0.12 


141 


10 09 00 


12 18 21 


18.09 


1.43 


7.91 


276.51 


3.97 


58.39 


65 


-1.76 


0.11 


143 


10 08 32 


12 18 34 


18.10 


1.32 


1.18 


292.99 


14.29 


14.10 


6 







Note. — This Table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 



Table 3 

Mean velocities and dispersions for different fields. 



Field 


(v) 
[ km s 


1] 


[kms 


cr 


GMOS, Field 1 


284.0 


± 


2.6 


7.4 


± 


3.4 


CMOS, Field 2 


283.6 


± 


2.8 


8.8 


± 


3.8 


GMOS, Field 3 


284.2 


± 


2.3 


12.5 


± 


3.6 


GMOS, all 


284.6 


± 


1.5 


9.9 


± 


2.1 


DEIMOS, all 


283.8 


± 


1.5 


9.9 


± 


2.1 


GMOS + DEIMOS 


284.2 


± 


1.0 


9.9 


± 


1.5 


GMOS + DEIMOS + M98 


287.0 


± 


0.9 


9.8 


± 


1.3 



Table 4 

Mass and central density estimates. 



Dispersion profile Sample "^^M^] [10« Mpfopc- 3] [{M/LyM 



Flat 8.7+i° 2.7 



+U.6 



2.6 ^-'-0.9 ^"-8 

Falling GMOS + DEIMOS S.ol^'^ 2.8l°-^ 24^^ 

Peaked 7.7^^^ 2.9lg| 23^^ 

Flat 7.9l^;* 2.8+^;^ 23+^ 

Falling GMOS + DEIMOS + M98 7-8l 2.8l°;^ 23^^ 

Peaked ^-^t^il 2.0lJ;^ 14^^ 
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Fig. 1. — Location of our targets centered on the Leol dSpli at (a, S) = (10:08:28, +12:18:18). Red giants observed with GMOS are shown 
as filled circles, DEIMOS targets are depicted as open squares. Overlaid ellipses designate respectively Leo I's core and nominal tidal radius 
at 3.'3, 12^6, respectively. The system's ellipticity and position angle are assumed to be 0.21 and 79°. 
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Fig. 2. — Radial velocity histograms of the GMOS data (solid line in the left panel), DEIMOS stars (dashed line, left panel) and combined 
GMOS plus DEIMOS set (right panel). All targets with velocities larger than 200kms~^ are shown. Within this range, two targets (at 251 
and 324 kms"'^, respectively) fall outside a ±3(T-cut. 
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Fig. 3. - Distributions of individual radial velocities (top panel) and mean radial velocities (bottom panels) as a function of their major 
axis distance. The symbols in the upper panel are the same as in Fig. 1, while the small dots additionally depict the 33 data points from 
M98. The solid horizontal line indicates the global mean of our sample, while the dashed lines denote the formal King radius at 12.'6 
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Fig. 4. — Tests for rotation in different subsamples of our data (as labelled): shown is the mean velocity difference with respect to rotation 
around axes at the respective position angles. A rotation signal around the major axis (at PA~ 79°) is distinguishable, but its significance is 
marginal. See text for a detailed discussion. 
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Fig. 5. — Velocity dispersion profile for different subsamples, as labelled, and rejection cuts. The numbers given refer to the sample size. 
Open circles refer to the full samples, whereas filled squares as well as the numbers in parentheses relate to the data cut at 3cr. 
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Fig. 6. — Top panel: Radial velocity dispersion estimator according to eq. (3), following the method of Walker et al. (2006b). Shown are the 
data for the combined GMOS + DEIMOS samples (small dots are the target locations), but the results do not alter after inclusion of the M98 
sample. The middle and bottom panels display contours of the statistical significances Phot a^nd Pcold fof the occurrence of any kinematically 
hot (middle) or cold (bottom) substructure. The individual panels show results for different choices of neighboring points included in the 
dispersion estimates. Contours are shown in intervals of 0.1. The dashed ellipse circumfers the nominal tidal radius and the two solid lines 
inscribe the interval of 3'-4', in which the radial dispersion profile exhibits a bump. See text for details. [See the electronic edition of the 
Jounal for a color version of this figure.] 
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Fig. 7. — Radial velocity dispersion profiles for the subsamples as in Fig. 5. From top to bottom each curve was calculated with increased 
binary fractions (labelled right of the outer bins). The solid curve and filled points refer to the observed profile assuming no binaries. Solid 
(open) squares and dashed lines are drawn from a maximization assuming an /(, of 0.4 and 0.8, respectively. In order to illustrate the eff'ect 
of a binary population on the errors on the velocity dispersion, the respective error bounds of the profiles with > are displayed as solid 
(open) triangles. 
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Fic;. 8. — Left panel: Surface brightness profile of Led from Irwin & Hatzidimitriou (1995) together with the fit of a Plummer law (solid 
line) and an exponential fit (dashed). The upper right panel shows fits to the dispersion profiles from the GMOS plus DEIMOS sample, 
and the analog after inelusion of M98's data is displayed in the lower right panel. Overplotted are fits of Plummer laws with different scale 
radii such as to reproduce a close-to flat (solid line) or falling (dashed) profile and also a modified profile, which incorporates the bump at 3' 
(dash-dotted). 
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Fig. 9. — Mass (upper panel) and density (lower panel) estimates from Jeans equation, assuming an isotrojjic velocity distribution. The 
left panels each refer to the GMOS + DEIMOS sample, whereas the right panels additionally account for the M98 data set. Different line 
types signify different approaches to fit the observed dispersion profile (see Fig. 8). The unphysical drop in the profiles of the left panels is 
due to the inclusion of an additional peak in the best-fit velocity dispersion profile, which exhibits a steep rising gradient in the inner regions 
opposed to a slower decline of the underlying Plumnier profiles. These counteracting gradients are not as distinct after inclusion of the M98 
sample so that the right panels do not display this drop. 
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Fig. 10. — Mass-to-light ratios for LG dSphs (after Wilkinson ct al. 2006a). The most recent individual mass estimates are from the sources 
cited in the text. The dashed line assumes a stellar (M/L)v of 1.5 in solar units and a dark matter component of SxlO^M©. Also shown for 
Leo I is the recent estimate of S06, who yield a lower value than the present work. 
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Fig. 11. — Predicted velocity dispersion profiles with varying degree of anisotropy for two different halo density profiles, namely, cuspy 
("NFW", top panel) and a cored halo (bottom panel). The inner regions are consistent with isotropic orbits, whereas the outer regions do not 
rule out a certain amount of radial anisotropy. The lower panel also contains a curve for the case of a radial varying anisotropy parameter, 
Pom, after Osiplfov (1979) and Merrit (1985), which represents our data reasonably well. 
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Fig. 12. — Mctallicity distribution for the 58 GMOS radial velocity member stars, for which EWs could be measured. These were convolved 
by observational uncertainties to yield the solid curve. Overplotted are a modified simple closed-box model (SCB, dashed line) and a model 
using Prompt Initial Enrichment (PIE, dash-dotted line), scaled to the same number of stars. 
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Fig. 13. — Spectra of the most metal poor star in our sample and a typical spectrum of a target around the peak metallicity. Apart from 
the three dominant Ca lines there are a few weak neutral metal lines detectable, which are not distinguishable from the continuum in the 
metal poorer target. 
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Fig. 14. — CaT metallicities of our radial velocity members vs. their elliptical radius (top). Nominal core and tidal radius are also denoted 
by vertical lines. The bottom left shows the cumulative number distribution versus radius of the metal richer component (solid line) and 
the metal poorer population (dotted line). These axe pra<;tically indistinguishable. MDFs at different radii are displayed in the bottom right 
panel. These density distributions were convolved by individual measurement errors. No apparent trend of [Fe/H] with radius is discerned. 



